Permanent confinement in compact QED 3 with fermionic matter 
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We argue that the compact three dimensional electrodynamics with massless relativistic fermions 
is always in the confined phase, in spite of the bare interaction between the magnetic monopoles 
being rendered logarithmic by fermions. The effect is caused by screening by other dipoles, which 
transforms the logarithmic back into the Coulomb interaction at large distances. Possible implica- 
tions for the chiral symmetry breaking for fermions are discussed, and the global phase diagram of 
the theory is proposed. 



It is well known that the pure compact electrodynamics 
in three dimensions (3D) is always in the confined phase, 
i. e. in the phase where magnetic monopoles form a 
neutral plasma of free magnetic charges yj. When the 
gauge field is coupled to matter, however, new possibili- 
ties emerge. In particular, Abelian gauge fields coupled 
to massless relativistic fermions (CQED3) arise naturally 
in several theories of high temperature superconductiv- 
ity 0- 0- CQED3 also represents maybe the simplest 
theory that should contain the physics of confinement 
and chiral symmetry breaking, and as such it may be used 
as a valuable toy model for the QCD. It has been argued 
that coupling to massless relativistic fermions turns the 
interaction between monopoles from 1/x to — ln(a;) at 
large distances x 0, [U, [g, 0, so that the deconfined 
phase with bound monopole-antimonopole pairs may be- 
come stable at low (effective) temperatures, in close anal- 
ogy to the Kosterlitz-Thouless (KT) transition. If so, 
this would suggest that, at least for some values of cou- 
pling constants, compactness of the gauge field may be 
neglected at large distances, and that the results pertain- 
ing to the continuum QED3 8] arc likely to remain valid. 
In particular, the chiral symmetry breaking instability of 
the continuum QED3 that has been linked to antifer- 
romagnetism in cuprates might then be expected to 
survive essentially intact. 

In this paper we consider the effects of finite density 
of monopole-antimonopole pairs (dipoles) on the puta- 
tive, KT-like, confincment-deconfincmcnt transition in 
the CQED3. We find that, in stark contrast to two di- 
mensions (2D), screening by logarithmically interacting 
magnetic charges bound to dipoles in 3D always alters the 
form of the interaction between the charges, from — ln(a;) 
at large distance, to 1/x. We demonstrate this first by 
an elementary electrostatic argument, and then support 
it by a controlled momentum shell renormalization group 
on the equivalent sine-Gordon theory. As a result, largest 
dipoles are always unbound, and would be expected to 
provide linear confinement of electric charges. Finally, 
we argue that the absence of the phase transition in the 
cQED 3 in terms of magnetic monopoles should imply 
the same in terms of fermions. Relying on the exist- 
ing numerical results and qualitative arguments, wc sug- 
gest that the phase of free monopoles corresponds to the 
(confined) phase of broken chiral symmetry for fermions, 
and that the (deconfined) chirally symmetric phase exists 



only in the continuum limit, and above a certain number 
of fermion species Q- The conjectured phase diagram is 
presented in Fig. ^ 

We are interested in the CQED3 defined on a cubic 
lattice: 

s = s f fcy-/j^co S (y, (i) 

where F^ v — e^peppa^-pOu is the lattice version of the 
electromagnetic tensor, and the coupling to fermions on 
a lattice may be written as 

N/2 

s f = 5 £ E fc {xye " {x) ^ {x + ^ (2) 

x,v 1=1 

using the standard staggered fermions. Here, x denotes 
the lattice sites and v the direction of the links. In the 
continuum, Sf corresponds to N species of four compo- 
nent Dirac fermions |lfj. 

The CQED3 has been previously studied numerically 
in fidj. Here, to make progress analytically and follow- 
ing [?]|, we approximate Q with 

^«^E[ F ^W-27rn^(x)](/3 + T ^) (3) 

[Fftuix) - 2irn liV (x)] , 

where 1/|A| should be understood as an inverse of the 
square root of the lattice gradient squared. The term 
proportional to the coupling f3 in Eq. @ is the standard 
Villain approximation to the second term in Q. The 
form of the term proportional to N in is motivated 
by the known continuum limit of the fermion polariza- 
tion: neglecting compactness of by setting n M „ (x) = 
would give the standard one-loop Gaussian action in the 
continuum, with the correct coefficient @- Just like in 
the ~ (3 term, here we also retained only the leading 
power of fields and their lattice derivatives. Eq. (3) may 
be understood as the quadratic approximation to CQED3 
in Eq. (1). It leads to the expected logarithmic interac- 
tion between the magnetic monopoles, and therefore ap- 
pears to contain the essential physics of CQED3 we wish 
to address. 
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FIG. 1: The proposed phase diagram of the CQED3. True 
phase transition between the chirally symmetric (CS) and chi- 
ral symmetry broken (CSB) phases of fermions occurs only in 
the continuum limit (f3 — » 00). Monopoles are in the plasma 
phase at all f3 < 00 and N < 00. Dashed line marks a 
crossover that corresponds to the N CTOBB (f3), as discussed in 
the text. 



Standard duality transformations 0] on © l ea d to 
the sine-Gordon action in the continuum 



(V$(x))|V|(V$(£)) 
tt 2 (N+ 16/3|V|) 



- 2ycos$(f) 



( 4 ) 

where y <C 1 is the fugacity of monopoles, and x is now 
the position vector in the continuum. In presence of 
fermions (N 7^ 0), at small momenta one can expand in 
the original gauge coupling /3 in Eq. (4), which becomes 
a coefficient of the term quartic in derivatives, and thus 
irrelevant by power counting. We may therefore safely set 
it to zero, and return to its (non-universal) effects on the 
phase diagram later. We also find it useful to have an al- 
ternative representation of the partition function defined 
by S s g, 



Xi exp 



i=l 



(5) 

in terms of unit magnetic charges qi = ±1 at an effec- 
tive temperature T = 4/N + 0{(3), interacting at large 
distance via V(x) = — ln(xA). A is the ultraviolet cutoff 
implicit in Eq. (4). The sine-Gordon field in Eq. (4) is 
then ($(x)) = {i/T)J2iqi{V(\S-Xi\)), where (...) dc- 
notes averaging with respect to Z. 

The excess of free energy due to a single isolated 



monopole in the sample of a linear size L is therefore 
1 



AF = -InL-TlnL 3 , 
2 



(6) 



so that for T < 1/6 monopoles would be expected to be 
bound in pairs. This corresponds to a "critical" number 
of fermion components N = 24 Q , , Q . This argu- 
ment neglects the effect of other dipoles, however. In the 
standard KT transition their presence does not change 
the critical temperature obtained by the simple energy- 
entropy argument, since the KT fixed point lies at zero 
fugacity [ijj. We show next that the situation in 3D is 
radically different. 

Let us assume a distribution of charges p(x) interacting 
via — ln(x) interactions in 3D, located in a region of a 
finite size R. At a distance x ^> R we can write the 
interaction as 



V(x) 



d 3 x'p(x') 



In \x\ 



(7) 



in spirit of multipole expansion [Llj . Thus, for a medium 
with a monopole density p{x) and a dipole moment den- 
sity P{x) — p{x)x , the potential is given by 



V(x) 



d 3 f (- In \x- x\)[p{x') - V' ■ P{x' ) + ...]. (8) 



The energy of a small dipole with a moment P = qr 
(q = 1) at x in a weak potential V(x) is still P ■ VV(x). 
Assuming non-interacting dilute dipoles at a temperature 
T then gives the average dipole moment at x 



- 2 ) 



(P($)) = (-W(f)). 

Since the density of dipoles is proportional to y 2 , 
polarizability of the medium is thus 



(9) 



the 



X 



3T y ' 



(10) 



where the constant of proportionality depends on the pre- 
cise form of the short distance regularization of the log- 
arithmic interaction. 

Writing Eq. (8) in Fourier space and combining Eqs. 
(9) and (10) into P{q) = — £x<zV(?), we finally obtain the 
interaction due to an external point charge Q in presence 
of a finite density of dipoles 



V(q) 



Q 



or, in real space 



V(x) 



1 1 



XT 



(ii) 



(12) 



This implies that in presence of other dipoles the first 
(energy) term in Eq. © scales like l/L, so that the en- 
tropy always dominates, and the largest dipoles should 
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lations (H} and JT^J is 



FIG. 2: The flow diagram in the y -a plane for N < N cr oss 
(a), and N > N cros3 (b). The non-trivial fixed point in (a) 
moved to lower right quadrant (not shown) in (b). 



dissociate into free monopoles. The 'dielectric' phase of 
logarithmically interacting monopoles in 3D thus, unlike 
in 2D, appears destabilized by dipole screening, as antic- 
ipated in jl5| . A similar conclusion has apparently also 
been drawn before in the study of instantons in quantum 
antiferromagnets [l6|. 

The renormalization group treatment of the equivalent 
sine-Gordon theory closely parallels the electrostatic con- 
sideration. Integrating out the modes of in Eq. Q 
with momenta A/6 < q < A one finds that the original 
couplings are renormalized as 



N(b) 

y(P) 

0(6) 



N + 0(y 3 ), 
6V~^ G>(0) +0(y 3 ), 

+ l y 2 e -G>(0) J d 3 x x 2 {e G > ( X )_ 1) 

-0(y 3 ), 



(13) 
(14) 

i 

(15) 



where a is the coefficient in front of the (VQ(x)) 2 term 
in S s g, Eq. (0J. Such a term is initially absent, but gets 
generated to the second order in y. Formally this is just 
what one expects, since there is no particular symmetry 
in the Eq. (@J that would prohibit its appearance. Physi- 
cally, as we have seen, this simply expresses the screening 
effect of smaller dipoles onto the interaction between the 
monopoles in a larger dipole. This term was effectively 
overlooked in the previous work on cQED 3 Q, @, 0, 
, where the authors consequently arrived at a different 
final result. 

The correlation function appearing in the recursion re- 



d 3 q 



n Ne iq ' x 



A/ b <\ ql <A (27T)3 2^3 +7T»W 



(16) 



and requires some care in evaluation. Assuming a low 
cutoff, or the effective temperatures TA 3 <C y, we may 
rewrite it as 



G>(x) 



7T 2 N 



2 A 



d 3 q 



(27T) 



q 2 



(A/6) 2 



n 2 Na 
„ 2 

f), 



(17) 



A 2 ' 



where the integration is now unconstrained, but we have 
introduced a smoothing function q 2 / (q 2 + A 2 ) instead of 
the sharp cutoff 01 • We have checked that when applied 
to the standard sine-Gordon problem in 2D, our regular- 
ization scheme gives the correct anomalous dimension at 
the KT transition. With this regularization we find 



G>{x) = 

and therefore 
N = 

y = 



1 TT 2 Nln(b) __ xA 
47r2 + 7r 2 Aa e 



0((ln&) 2 ), 



+ O(f), 



ttN 



8(2 + 7r 2 Aa) 
8tt 2 A 



2 + TT 2 Na 



0(f), 



(18) 

(19) 
(20) 

(21) 



where y — y/A 3 and a = a/ A, are the dimensionless 
couplings, and z = dz/dhxb\b=i- 

Two important features of the flow equations 119fl - 121|) 
should be noted: 1) the relevant coupling a, even if ab- 
sent initially, becomes generated at lower cutoffs, and 2) 
N is marginal to 0(y 3 ). We suspect N to be an ex- 
actly marginal coupling, since the action for slow modes 
should be analytic in (q/A) 2 , so the coefficient of the 
non- analytic \q\ 3 term can not get renormalized to any 
order |l8|. 

The flow in the y 2 -a plane for N < N CTOSS and N > 
Across, A cross = 48/-7T, is depicted in Fig. [3 Starting from 
a = 0, y = yo, the fugacity in the former case monoton- 
ically increases as b — > oo, while in the latter it begins 
to increase only for b > b* , with b*(y ) w l/(967ry ) at 
N ~S> 1. Nevertheless, since y = is an invariant line un- 
der the scale transformations, fugacity must increase at 
large enough b for any initial value. We interpret this up- 
ward flow of the fugacity as an indication that monopoles 
are always in the plasma phase. 

Interestingly, massless fermions still make themselves 
felt in such a plasma, since the non-analytic |V| 3 term 
in Eq. (4) translates into a power-law behavior of the 
screened interaction at long distances. Instead of the ex- 
pected Debye-Hiickel exponential decay, in the quadratic 
approximation to (4) we find Vs C r(a;) ~ — Q/(ir 4 y 2 N\x\ 6 ) 
at large x. The reader may recall a similar phenomenon 
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of Friedel oscillations in metals, where a non-analyticity 
due to the Fermi surface produces a power-law behav- 
ior of the screened interaction as well. In the present 
case the sign of the above power-law term indicates over- 
screening, presumably due to the extremely long range 
of the bare (logarithmic) interaction between monopoles. 

Let us next turn to possible implications of our find- 
ings for fermions. First, the dimensionless coupling 
(3=1/ (e 2 l), where e is the dimensionfull charge, and I the 
lattice spacing. The lattice gauge field O^x) = elA /1 (x), 
where A^(x) is the gauge field in the continuum. The 
continuum limit is achieved by keeping e constant and 
taking / — > 0. In other words, (3 — > oo should naively 
correspond to the continuum QED3. In this theory chi- 
ral symmetry for fermions is expected to be broken, and 
the condensate (jffi) ^0; for N < N c h, where N c h is of 
order unity @, E3i E3- On the other hand, for (3 < 00 
we argued that once fermions are integrated out, the the- 
ory has no phase transition. If correct, this implies that, 
had we proceeded in reverse and integrated out every- 
thing but the fermions, we would have had to find them 
in a single phase at any finite (3. Relying on numeri- 
cal evidence that free monopoles seem to enhance chiral 
symmetry breaking we propose that it is the phase 
with broken chiral symmetry that survives a finite (3. 

It was found previously |2l| that the critical N for chi- 



ral symmetry breaking at j3 — is significantly larger 
than iV^ in the continuum limit. We expect that for 
(3 = 0, there is no true phase transition in the CQED3, 
just the crossover occurring at N CTOSS . Making the sim- 
plest assumption that (3 always stays irrelevant, the same 
would hold at any (3 < 00, but with the crossover line 
shifted to iV cross (/3) = Across - 16/3A + 0{{(3kf). The 
conjectured phase diagram that summarizes this discus- 
sion is presented in Fig.^ (^^) 7^ everywhere, except 
on the line 1/(3 = 0, N > N ch . 

To conclude, we argued that magnetic monopoles 
are always free in the CQED3 with gapless relativistic 
fermions, and the theory is thus expected to be perma- 
nently confining. We proposed the phase diagram, and 
conjectured that chiral symmetry for fermions is always 
spontaneously broken. Our results seem to be in quali- 
tative agreement with the early numerical results of ref. 
[2^|. They can also be generalized to the case with non- 
rclativistic fermions with Fermi surface [23j . Although we 
have studied only the U(l) theory coupled to fermions in 
this paper, similar arguments should be applicable to the 
case of critical scalar field Q , ■ 
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